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Abstract. In a series of papers llLombardi fc Schn eider 200l], I2002T) we studied in detail the statistical properties 
of an interpolation technique widely used in astronomy. In particular, we considered the average interpolated map 
and its covariance under the hypotheses that the map is obtained by smoothing unbiased measurements of an 
unknown held, and that the measurements are uniformly distributed on the sky. In this paper we generalize the 
results obtained to the case of observations carried out only on a finite field and distributed on the field with a 
non-uniform density. These generalizations, which are required in many astronomically relevant cases, still allow 
an exact, analytical solution of the problem. We also consider a number of properties of the interpolated map, 
and provide asymptotic expressions for the average map and the two-point correlation function which are valid 
at high densities. 
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1. Introduction 

Interpolation techniques play a central role in many phys- 
ical sciences. In fact, experimental data can often only 
be obtained at discrete points, while quantitative, global 
analyses can normally be performed only on a field. A 
classical example of such a situation is given by meteoro- 
logical data, such as temperature, pressure, or humidity: 
these data are collected by a large number of ground-based 
weather stations, and then need to be interpolated in or- 
der to obtain a continuous held. 

The situation is, apparently, very similar in 
Astronomy. Indeed, many astronomical observations 
are carried out "discretely," i.e. data are available only 
on some locations of the sky (typically corresponding to 
some astronomically significant objects, such as stars, 
galaxies, quasars). If there is some reason to think that 
the data represent discrete measurements of a continuous 
held, then the observer will want to interpolate the data 
in order to obtain a smooth map of the quantity being 
investigated. 

In reality, astronomical observations have a character- 
istic that make them quite peculiar with respect to other 
physical experiments: in most cases, it is not possible to 
choose where to perform the measurements. A meteorol- 
ogist, for example, can always decide to put his weather 
station in a particular location, or to have weather sta- 
tions regularly spaced; this, clearly, is impossible for an 
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astronomer. As a result, it is sensible to perform a sta- 
tistical analysis of interpolation techniques by consider- 
ing the measurement locations as random variables, i.e. 
by performing an ensemble average on the positions of 
the astronomical objects used for the analysis (this tech- 
nique has been already used by several autho rs; see, e.g., 

lLombardi fc BertmlllQQSUvan Waerbekell2000lh 

In a series of previous p apers (|Loanb^Tdj_^L-Scii neid^ 

l200ll hereafter Paper I, and lLombardi fc Schneiderll2002L 
hereafter Paper II), we analyzed the statistical properties 
of a widely used interpolation technique. In particular, we 
considered a set of measurements {f n } performed at loca- 
tions {On}. The measurements were taken to be unbiased 
estimates of a field f(9) at the relative positions, i.e. 



{fn) = /(*«) , 



(1) 



where the brackets (•) denote the expectation value of the 
enclosed quantity. The discrete measurements {f n } were 
then interpolated using the following technique. First, we 
introduced a positive function w(cf>), which describes the 
"influence" of measurements performed at 9' = 9 + <fi on 
the interpolated field f{9). This field was defined as 



n ) 



J2tM0-0n) 



(2) 



where N is the total number of observations. Equation (J5J), 
indeed, is a standard i nterpolation technique (see, e.g. 
ICressielll993l lLanJll98.il called "moving weights," "mov- 
ing average," or "distance weighted average" [this last 
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name is due to the fact that normally the weight 
function w used in Eq. (0) depends only on the dis- 
tance ||0 — 9 n \\). Although other techniques are clearly 
available (see, e.g. iBernardeau Sz van de WevgaertJllflQfit 
ISchaap k van de Wevgaerdl200ri lLombardill2002(r this is 
probably the interpolation method most often used in 
Astronomy. 

In this paper we study the expectation value and the 
two-point correlation of the smoothed map f(9) under the 
hypotheses that: 

1. The measurements {f n } are unbiased estimates of the 
field / [cf. Eq. QJ] with errors. The errors e n = f n — 
f(9 n ) ar G taken to be independent random variables 
with vanishing mean (this is clearly equivalent to the 
unbiasness of {/«.}) : 



<e„) = , 



(3) 



2. The measurement locations are independent random 
variables distributed according to a know density field 
p{9) inside a given observation area Q (i.e., are a non- 
homogeneous Poisson process). 

Hence, this work generalizes the results obtained in 
Paper I and II in three directions: (i) a non-constant vari- 
ance a 2 (9) is considered; (ii) the density of measurement 
locations p can vary on the field; (m) no restriction is 
put on the size of the observation area f2, which can be 
finite. Surprisingly, although these generalizations signifi- 
cantly widens the applicability of our results, they do not 
make our method more complicated at all. Indeed, as we 
will see below, the problem seems to find a very natural 
description in the more general framework used here. 

It should be stressed that the generalizations carried 
out here are very important in the astrophysical context. 
Indeed, astronomical data will normally be available only 
on a limited area of the sky, and so boundary effects have 
to be taken into account. Moreover, often the measure- 
ments will not be uniformly distributed on the observed 
area. This happens, for example, for data based on stars, 
which have a higher density when one approaches the 
galactic equator. However, even for astronomical objects 
that arc, in principle, uniformly distributed on the sky 
(e.g., distant galaxies or quasars), we might need to deal 
with a non-uniform distribution because of observational 
effects (e.g., because of a non-constant sensitivity on the 
field of the detector used, of dithering patterns, or of the 
presence on the field of bright objects that interfere with 
the measurements). 

The paper is organized as follow. In Sect. 2 we carry 
out the various generalizations in turn. The properties of 
the average and the two-point correlation function of the 
smoothed field are discussed in Sect. 3. Finally, we sum- 
marize the results of this paper in Sect. 4. 



2. Main results 

2.1. Position dependent weight function 

Looking at Eq. we can note that the interpolation 
point 9 does not directly enter the problem, but is basi- 
cally only used to "label" the interpolated point. For our 
analysis, indeed, it is convenient to fix that point to, say, 
9 = 9a, and to rewrite Eq. (J2J) as 



_ J2n=lfnWA(9 n ) 



(4) 



where we called /a = J{9a) and wa{9') = w(9a — 9'). 
In this new notation, in the hypotheses used for Paper I 
(homogeneous Poisson process with uniform density p 
for measurement locations, infinite field), the expectation 
value {f a) can be evaluated exactly from the equations 



Q A (s)=p [ [ e -" A W-l] A 2 9, 
Jn 



Ya{s) = exp[Q j4 (s)] , 
1 



Ca(wa 



s Y A (s)ds , 



(U) = P / f(9)w A (9)C A (w A (9)) d 2 e , 
Jn 



(5) 
(6) 
(7) 

(8) 



where fl is the observation area (taken to be very large 
compared with the typical scale of the weight function 
wa) and Pa is the probability of having no point inside 
the support of wa (i-e. Pa = exp(— pita), with tta the area 
of the support of wa)- 

Equations show explicitly that the location 9 

does not enter our problem. As a result, looking again at 
Eq. ©, there is no need to take the weight function to 
be of the form w{9 — 9 n ), but we can instead consider 
in the same framework the more general form w(9,9 n ). 
As we will see below, the trivial generalization described 
in this subsection is actually a fundamental step for more 
interesting results. 

2.2. Finite fields 

We now focus on a slightly different generalization, namely 
the use of finite fields. We observe that having no data 
outside a given field is totally equivalent to having data 
everywhere and using a vanishing weight for locations out- 
side the field. In other words, even if our observations are 
confined on a small part of the sky, we can always imag- 
ine to have data in the whole sky, by generating arbitrary 
values for the data locations and values, and then discard 
the arbitrary data by using a vanishing weight for them. 
In turn, from the form of the integrands in Eqs. © and 
© , we see that the integrals actually need to be performed 
only on the domain of the weight function (the integrands, 
indeed, vanish ot the points where wa vanish). As a re- 
sult, Eqs. (JSHEJ) can still be used, provided we interpret SI 
as the observation field. 
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2.3. Non-uniform density 

We can now finally generalize Eqs. JHHE1) to non-constant 
densities. We first observe that the remaining spatial vari- 
able 8' that appears in the definition J3J is a dummy 
variable. Indeed, 8' is basically used to "name" locations, 
but does not really play any role on interpolation process. 
For example, performing an arbitrary bijective (i.e., one- 
to-one) mapping 8 i— > rj described by the function rj(8), 
will not change the value of /a, provided that we use the 
"mapped" weight function w A \6) = w^(ji{8)). On the 
other hand, when performing the transformation 8 t— > rj, 
we are bound to change the density distribution of objects. 
More precisely, if the objects are uniformly distributed on 
the plane rj with density p^) ; they will be distributed ac- 
cording to a non-uniform density p^ (8) on the 8 plane. 
The final density, indeed, is given by 

fdri(d)\ 



(0) = 



det 



V 96 J 



(9) 



This observation suggests a possible way to study non- 
uniform densities. Suppose that we intend to study the 
expectation value of Ja of Eq. Q when the locations on 
the 8 plane are distributed according to a non-uniform 
density p^{8). Then, we can look for a one-to-one map- 
ping 6 i— ► rj, such that the corresponding density p^ v \ 
evaluated according to Eq. @, is uniform; we also mod- 
ify the weight function accordingly. Now, since in the r\ 
plane the locations are uniformly distributed, we are in 
the position of studying the problem using the technique 
developed in Paper I. Moreover, since the value of Ja does 
not depend on the coordinates r\ or 8 used, we can directly 
use the result obtained. 

The method described above clearly allows us to solve 
a much more general problem but it has also two main 
problems. From the theoretical side, one has to show that 
it is possible to find a one-to-one mapping that satisfies 
our needs (namely, that p^ is uniform) . From the practi- 
cal side, it might be non-trivial to find the function r)(6); 
moreover, for every point 8 a, one needs to transform the 
weight function w A (8) = w{8 a, 8) into a weight function 
on the r] plane (see above). We now address both prob- 
lems, showing that our equations can be reformulated in 
a way that naturally allows for non-uniform densities. 

First, we explicitly show that, for every density distri- 
bution p( e \d) it is always possible to find a one-to-one 
function rj{6) such that the corresponding p^, evaluated 
from Eq. (0, is constant. Let us, in fact, consider the 
transformation 



rim 



(0) 



9' 1 ,e 2> ... ) e M )de' 1 if * = 1 , 

otherwise , 



(10) 



where M is the number of dimensions of 8 (typically, in 
Astrophysics, M = 2). This transformation satisfies the 
following properties: 

— It is continuous provided p^ is continuous. Situations 
where p^ is non-continuous will be considered at the 
end of this subsection. 



It satisfies by construction Eq. ©. 
Jacobian matrix drj/dd has the form 



dd 



1 
1 



*\ 






V 0... 1/ 



Indeed, the 



(11) 



where the stars (*) denotes uninteresting terms. 
— If // e ) is strictly positive, the transformation is one- 
to-one. Indeed, if p^ > 0, then r](8) is injective and 
also (clearly) subjective if restricted to its codomain. 
In case p^ vanishes on some points, we can always 
use the argument of Sect. 12.21 In other words, a van- 
ishing density in a region implies that the probability 
of having points in this region vanishes as well; on the 
other hand, having no measurements on that region in- 
troduces a finite field described in Sect. 12.21 Hence, we 
can assume that, inside f2, the density is strictly 
positive. 

This proofs the existence of a one-to-one function r){6) 
with the required properties (note that there are many 
possible choices for rj(8), but they are all equivalent for 
our purposes). 

We now turn to the second problem, namely the prac- 
tical difficulties in applying the technique discussed in this 
section. Suppose again that we are interested in evaluating 
the expectation value of j 'a of Eq. (0J with a non-uniform 
density p^ e \d). Then, we can use Eq. I|10|) to convert the 
problem into the rj plane, so that the corresponding den- 
sity is unity. We can then finally apply Eqs. (jHHEl 011 V, 
using p^ 1 ' = 1. In particular, for Eq. © we have 



Qa(s)=p m 



■n(Q) 



1 



d 2 v 



e-^W-ll p^(6)d 2 9. 



(12) 



Note that in the second line we have operated a change of 
variable from r\ back to 6; note also that we have used the 
fact that w( 6 \6) = w^(r)(6)). Hence, we can still use an 
equation very close to Eq. (JSJ: we just need to integrate 
also over the variable density p( \ Similarly, for Eq. JSJl 
we find 



{Ia)=P M 



7,(0) 



f^(r,)w^( V )C A (w^(r,)) d 2 r, 



We 



f^(6)wf(6)C A (wf(6)) 

= / f { A ) {e)wf{d)c A { w ^\6))p^{e)d 2 e. (13) 

Equations (JSJ and (JJJ remain unchanged. Hence, by sim- 
ply using a transformation of variable back to 6, we have 
been able to obtain a solution of the problem that does 
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not involve rj any more. This shows once more that the 
problem, as expected, does not depend on the details of 
the choice of the transformation 9 77 (see comment at 
the end of the previous paragraph). Hence, in the follow- 
ing we will drop the superscript {9) used in this section, 
and hence we will always assume that all functions are 
evaluated in the coordinate system defined by 9. Finally, 
note that the expressions in the last lines of Eq. 1(12(1 and 
HI 3(1 are still valid for vanishing densities; in other words, 
we can either include the finite-field effects on the defini- 
tion of f2, or just put p^ e \0) = outside the observation 
field. 

Before closing this subsection, we consider the case of 
a non-continuous density field p^ e \ We recall, indeed, that 
with Eq. 1 (10( 1 we have been able to provide a one-to-one, 
continuous transformation rj(9) only under the hypothesis 
that p^ be continuous. Although this hypothesis is not 
needed, it is non-trivial to exhibit a mapping with the 
require properties in the general case of a non-continuous 
density. In reality, this problem is only apparent. Note, in 
fact, that the density p^ enters Eqs. ((13(1 and 1(15(1 only 
as a term inside an integral and thus the continuity of 
this function does not play any role in our problem. For 
example, if we convolve a discontinuous density with 
a Gaussian, 



i e Ho) 



1 



2nd- 



■ exp 



\e 



M12 



2a 2 



p( e \0-0')d 2 6' , (14) 



(0) 

we obtain a smooth function ps (9). This function, then, 

can be used in Eqs. 1(131) and ((15(1 at the place of the den- 
(0) 

sity. Note that, since ps (9) is smooth for any a > 0, we 
can apply the transformation 1(10(1 without any problem. 
Finally, we take the limit a — ► + , so that the results of 
the integrations l(13|) and 1(15(1 are not modified by the use 
of pi°^ instead of p^ . 




— 1 — 

WAP 
Pi 
pi 
P3 



Fig. 1. The effect of a non-constant density on the effec- 
tive weight. The plot shows, in the ID case, the quan- 
tity p(x)w c g(x) for three different densities, pi(x) = 1, 
P2(x) = 1 + H(x)/4, and p${x) — 1 — (cosx)/2. In all cases 
the original weight function has been chosen to be such 
that the combination w(x)p(x) is a Gaussian (see solid 
line plot). 
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2.4. Average map 

We summarize here the results obtained in this section. 
We have shown that the expectation value of /a can be 
evaluated from the set of equations 



Q A (s)= [ [e- SWA ^ ~l]p{9)A 2 6 , 



Ya(s) = exp[Q A (s)] , 
1 



Ca(wa 



I -P. 



A Jo 



'Ya(s) ds 



(U)= / f(9)w A (9)C A (w A (9))p(9)d 2 , 
where the probability Pa can be evaluated from 



Pa = exp 



p{9)d 2 e 



(15) 
(16) 
(17) 

(18) 



(19) 



Appendix A reports an alternative method to derive this 
set of equations. In the rest of the paper we will assume, 



Fig. 2. Effective weight function in presence of bound- 
aries. Three Gaussian weight function (shown in solid 
lines) centered on different part of the field £1 = [—5,5] 
produce significantly different effective weights (dashed 
lines). The weight function have been normalized accord- 
ing to Eq. J30I); p(x) — 1.5 is constant for this plot. 



without loss of generality, that the weight function satis- 
fies, for each point 9, the (position dependent) normaliza- 
tion property 



w{9, 9')p{9') d 2 6' = / w A (9')p(9 , )d 2 8' = 1 



(20) 



Indeed, since only relative values of the weight function are 
important in our problem, we can always suppose to deal 
with weight functions normalized according to Eq. 120(1 . 

Similarly to what was done in Paper I, we call w c sa — 
waCa(wa) the effective weight function, so that we can 
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write Eq. (JTSJl as 

(Ia) = I f(6)w cHA (0) P (0)d 2 6 
Jn 



Tp3 = (/a>(/s) 



(21) 



Note that, in contrast to Paper I, in the definition of the 
effective weight function we have not included the density 
p(6), which thus must be explicitly added in the integra- 
tion of Eq. J2J. This is convenient, because this way the 
effective weight is a function of the value of the weight 
function and not (directly) of the position 9, and because 
this way, as we will see in Sect. 13.11 the normalization 
property of the effective weight function is similar to the 
normalization (|20() of the original weight. The effects of a 
non-constant density and of a finite fields on the effective 
weight are shown in Fig. and [3 

2.5. Final solution: two-point correlation function 

We now turn to the generalization of the results of Paper II 
concerning the covariance of /, i.e. its two-point correla- 
tion function. Since the generalization procedure closely 
follows the one used in Sect. [5] for the average, we skip 
here many details and mainly report the final result only. 

Wc first recall that in Paper II we have defined f B and 
w B similarly with Ja and wa, with the only difference 
that now these quantities are calculated with respect to 
a different point 9 B . We then have defined the two-point 
correlation function of / as 



Cav(f;d A ,9 B ) = {}aIb) - (h)(h) , 



(22) 



and we have shown that this quantity is composed of two 
terms, Cov(/; 9a, 8b) = T a + Tp, where T„ is the noise 
due to measurement errors and Tp = Tpi + 7p2 — Tp3 is 
the Poisson noise [split, in turn, of three terms; see below 
Eqs. 

Using an argument similar to the one adopted in 
Sect. |2 we can generalize the results of Paper II to the 
hypotheses discussed in the items of Sect. ^ We show 
here only the final results and skip the proof, which is a 
trivial repetition of what was done above for the average. 

Q(s A ,s B )= f [e- SAWA ^- SBWB ^ -l]p(e)d 2 e , (23) 
Jn 

Y(s A , s B ) = exp[Q(s A , s B j] ■ (24) 

pOC poo 

C{w A ,w B )=v ds A ds B e- s ^- SBWB Y{s A ,s B ) , 
Jo Jo 

(25) 



T a = [ d 2 6p(e)<j 2 (9)w A (0)w B (e)C(wA(e),w B (0)) , 



(26) 



Tpi= / d 2 0p(9)[f(9)] 2 w A (e)w B (9)C(w A (9),w B (9)) , 
Jn 



(27) 



T P2 = [ d 2 e lP {9 l ) f d 2 6 2 p(9 2 )f(9 1 )f(9 2 )w A (9 1 )w B (9 2 ) 
Jn Jn 

x C(w A (9 1 ) + w A {9 2 ),w B (9 1 ) + w B {9 2 )) , (28) 



d 2 d x p(9 1 )f(9 1 )w A (9 1 )C A (w A (9 1 )) 



d 2 9 2 p(9 2 )f(9 2 )w B (9 2 )C B (w B (9 2 )) 



(29) 



We have called in Eq. v = 1/(1 - P A - P B + P AB ), 
where the probabilities Pa and P B can be evaluated from 
Eq. (JT5J, and Pas (the probability of having no points 
inside both the supports of wa and w B ) is given by 

P AB = cxp - / p{9) d 2 . (30) 



3. Properties 

In this section we will consider some interesting properties 
of the average (Sect. I3~T|) and of the two-point correlation 
function (Sect. l3~2|) of f(9). Hence, here we basically gen- 
eralize Sect. 5 of Paper I and Sect. 6 of Paper II. 

3.1. Average map 
3.1.1. Normalization 

By construction, for a constant field f(9) = 1 the 
smoothed function f A defined in Eq. returns on av- 
erage 1, a property related to the normalization of the 
effective weig ht rsee lLornharcTill20n2h . Indeed, if f{9) = 1, 
we find 



w eBA (9)p(9)d 2 9 



w A (9)C A (w A (9))p(9)d 2 e 



= - X p / dse QA{s) f d 2 9 p{9)w A {9)e~ SWA{e) , 
1 — "a Jo Jn 

(31) 

where Eqs. (|15H18|) have been used in the last step. The 
last integral is just —Q' A (s), so that 



1 



1-PaJo 

1 n 



Q' A (s)e QA{s) ds 



I -Pa 



a(s) 



(32) 



Hence, the effective weight has the same normalization 
property as the weight function wa [see Eq. (|2U|l ]. 

3.1.2. Scaling 

Suppose we rescale the weight function w A (9) into 
k 2 w A (k9), and at the same time the density p(9) into 
k 2 p(9); then we can verify using Eqs. \ 1511 18(1 that the ef- 
fective weight is rescaled similarly to wa, i-e. w c sa{9) i— > 
k 2 w cS A(k9). 
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This scaling property suggests the following definition: 



M A 



{I -Pa) [ [w A (0)] 2 p(0)d 2 
Jn 



(33) 



This quantity represents the number of "relevant" points 
used in the average, i.e. the expected number of locations 
for which the weight w A is significantly different from zero. 
The (1 — Pa) term in Eq. I|33|l is introduced in order to 
compensate for finite-fields effects [cf. the similar factor in 
Eq. (J7J)]; for example, for the top-hat function, this guar- 
antees that Ma — pnw A / [l — exp( ~P^w A )] > 1 for any 
density. In any case, the above definition finds its main 
justification from the properties that the quantity Ma so 
defined enjoys (see below). Following Paper I, we call Ma 
the weight number of objects; similarly, we define the ef- 
fective weight number of objects M cBA using the effective 
weight w c sa instead of wa in Eq. (|33l) . Note that the re- 
lated weight area A and effective weight area A c g, also 
introduced in Paper I [see Eq. (39) there], cannot be de- 
fined for the general case of a non-uniform density. 

3.1.3. Behavior of waCa{wa) 

A study of w A C A (w A ) can be carried out using the same 
technique adopted in Paper I. Since Ya(s) > for every s, 
C A (w A ) decreases as wa increases. Regarding wa Ca(wa), 
from the properties of Laplace transform (see, e.g.. lArfkenl 
Il985l see also Appendix D of Paper 11) we have 



WcBA = w a C a {wa) = Y A (0) + C[Y A ](w A ) , 



where £[•] indicates the Laplace transform. Since Y A (s) = 
Q' a (s)Ya(s) < 0, wc find that w c g A (w A ) increases with 
wa- This implies that there must be a value w A of the 
weight function wa such that C A (w A ) > 1 if wa < wa, 
and Ca(wa) < 1 if wa > wa- Indeed, since Ca is mono- 
tonic, the equation w e ff A (w A ) = wa can have at most one 
solution; however, this equation must have at least one 
solution because both wa and w e s A satisfy the same nor- 
malization property [cf. Eqs. l{20*|l and 432)]- The quantity 



D= [w A (0) + w cfiA {0) - 2w A ] 
Jn 



[w A (0)-w cBA (0)]p(0)d 2 8>0 



(35) 



is positive or null because the integrand is non-negative 
everywhere. By expanding the integrand we find 

0<D= [ [w A (e)] 2 p(0)d 2 0- [ [w cBA (8)] 2 p(8)d 2 6 
Jn Jn 

-2w A f [w A {e)-w eSA (0)]p(e)d 2 e 

Jn 



m 

1 ( 1 



1 



1 - P A \Ma M c 



cir.4 



where the normalization of w A and of w e s A has been used. 
Hence, we find M e SA > Ma- 



We now consider the limits of w c s A (w A ) for small and 
large values of w A , 



Y(s) I 
lim w A C A (w A ) = lim — = — 

wa^oo s ^o+ 1 — Pa 1 - Pa 



(37) 



Since w c sa(wa) is monotonic, (I — P A ) is a superior 
limit for the effective weight function. This property, in 
turn, can be used inside the definition of M cBA to obtain 



MJ A = (1-P A ) I [w A (0)] 2 p(6)d 2 9 
Jn 

< f w A (0) P (e)d 2 e = i. 

Jn 



(38) 



In other words, the effective number of objects will always 
exceed unity, independently of the weight function used. 
Regarding the other limit we have 

lim w A C A (w A ) = lim - — ] ——Y A (s) = - Pa . (39) 

W A ^0+ S^QC 1 — t' A 1 — P A 

Hence, since w b sa(wa) vanishes if wa = 0, the effective 
weight has a discontinuity at if P A ^ 0. 



3.1.4. Limit of high and low densities 

At high densities (p — > oo) only values of Qa{s) close to 
s = are important, because for large s, Ya(s) vanishes. 
Hence, we expand Qa(s) by writing 



(34) Q A (s)=J2 



fe=l 



(-l ) k s k S Ak 
k\ 



where S A k is the kth moment of w A : 

s Ak ee f [ WA (e)] k p(9)d 2 d . 

Jn 



(40) 



(41) 



The normalization (|20|l implies S A \ — 1, and so to first 
order Y A (s) ~ e~ s . We have then 



C A (w A ) 



1-PaJ 



e - SWAe S ds 



1 1 

l-P A l+w A 



(42) 



In the limit of low densities (p — > + ), instead, 
Y A (s) —> 1 and 



C A {w A ) 



1 1 



(43) 



1 Pa w a 

Expanding Eq. i|19|) to first order in p we find, for wa > 



(36) w AcS = w A C A {w A ) 



p(8)d 2 d 



(44) 



Hence, the effective weight converges to a top-hat function 
normalized to unity. 
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Fig. 3. The variance T a evaluated at different positions 
xa = xb- For this plot, we used a Gaussian weight func- 
tion with variance 1/2; the measurement error a was 
kept constant on the field. Note that, in agreement with 
Eq. (jSTJ, the quantity T a /a 2 is always larger than l/A/" e ff; 
curiously, the quantity 1/7V gives a good first-order ap- 
proximation for T/y even in the complex situation shown 
here. 

3.1.5. Moments expansion 

At large densities p, we can expand Ca{wa) in terms of the 
moments of wa defined in Eq. I|41|l . Calculations are basi- 
cally identical to the one provided in Paper I [see Eq. (66) 
of that paper] , with only minor corrections due to the dif- 
ferent definition of Ca- Hence, we skip the derivation and 
report here only the final result (up to the fifth term): 



(1 - P a )Ca(w a ) 



1 



SA2 



Sa3 



l + WA 
Sa4 ~ 



(l + w A ) 3 (1 + wa) 4 



+ 



(45) 



(1 + wa) 5 

3.2. Two-point correlation function 

The generalization of the properties of the covariance 
terms T a and Tp is, in most cases, trivial and closely fol- 
lows the generalization carried out in the Sect. l3.il Hence, 
here we will skip much of the proofs and just outline the 
main results. 

3.2.1. Normalization 

It can be shown that the Poisson noise satisfies a sim- 
ple normalization property: if f{9) = 1 is constant, then 
Tpi + Tp2 = Tp3 = 1, and thus the Poisson noise vanishes. 
A proof of this property can be carried out either with 
the technique described in Paper II, or, more easily, u sing 
the following argument, taken from lLombardil (2002). If 
f(9) — 1 is constant on the field, we will on average mea- 



sure (/„) = 1 for each point. Let us now assume for a mo- 
ment that we do not have any measurement error, so that 
/„ = 1 for every n. Then, we will always measure f{9) = 1, 
and thus (f(9 A )) = (f(0 B )) = (f(9 A )f(9 B )) = 1. In 
this case, thus, we find Tpi + Tp2 = Tp$ = 1. The situ- 
ation is actually the same even if the measurements are 
affected by errors: These, in fact, appear only in the eval- 
uation of To-, and thus the Poisson noise is left unaffected. 

3.2.2. Small and large separations 

In the limit 9a = the expressions for T a and Tp take 
a particularly simple form. Indeed, since wa = wb, to 
evaluate T a , Tpi, and Tp2 we just need C(wa,wa)', this 
quantity, in turn, can be easily shown to be C(wa, wa) = 
—C' A {w A ), where Ca (wa) is given by Eq. l(T7|) . 

If instead \9a — 9 B \ is large compared to the scale 
lengths of the weight functions wa and wb , then 



Q(sa, sb) 
Y(s A , s B ) 
C(w A ,w B ) 



Qa(sa) + Qb(sb) , 
Ya(sa)Y b (sb) , 
Ca(wa)C b (w b ) ■ 



(46) 
(47) 
(48) 



The following argument shows that in general 
C(wa,wb) > Ca{wa)Cb{wb)- First, observe that, since 
Pab > PaPb, one has v > (1 - Pa) _1 (1 - Pb)" 1 ■ 
Moreover, it can be shown that Q(sa,sb) > 
Qa{sa) + Qb(sb)'- indeed 



Q{s A , s B ) - Qa{sa) - Qb{sb) 

-s AWA (G) _ jj ^ c -s B w B (e) _ j-j^g) d 2 Q > q 



= /[e 
Jn 



(49) 



Finally then 

C(w A ,w B ) - Ca(w a )Cb{w b ) 

1 1 f°° Z -00 

t ±A t ±B Jo Jo 



x e 



^Q(s A ,s 



i(sa) + Qb(s B ) > Q 



(50) 



where the last inequality is a consequence of the mono- 
tonicity of the exponential function and of Eq. 149f) . 

3.2.3. Behaviour of T CT 

The normalization of the Poisson noise terms derived in 
Sect. 13.2.11 can be used to derive an upper limit for T a . 
Indeed, from the expression of T„ one sees that this quan- 
tity is very similar to Tpi, provided we replace f(9) with 
cr(9). On the other hand, from Sect. l3.2H1 we know that, if 
f(9) = 1, then Tpi < 1, because in this case Tpi+Tp2 = 1 
and Tp2 is positive. Hence we find T a < cr^ lax , where cr^j ax 
is an upper limit for <r 2 {9). 
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A lower limit for T a can be obtained from the inequal- 
ity J5DJ: 



T a > / a 2 (0)w A (9)w B (e)C A (w A (9))C B (w B (e))p(9)d 2 



(51) 



= f a 2 {6)w cSA (o)w cS B{e)p{e)& 2 e 

Jn 



Note that if 9 A = Ob, we have w c ff A (9) = w e SB{0) and 
the last integral in Eq. I|51|l is closely related to N e SA = 
•A/'offs- The two limits on T a discussed in this section are 
exemplified in Fig. [21 

3.2.4. Limit of high and low densities 

At high densities we can expand Q(s a ,sb) as done in 
Sect. EH for Q A : 



s A ,s B ) = J ^ -^-j[-s A w A (9)Y [-s b wb{0)] J - 1 



p(9)d 2 



-S A - S B 



(52) 



where in the last step we have retained only the first terms 
of the sum, and used the normalization of w A and wb [cf. 
Eq. (|2"0j) ]. Hence, in this limit 



C(w A ,w B ) ^ v 



1 



1 



1 + w A 1 + wb 



(53) 



Note that, because of the normalization (|2U[1 . both weight 
functions w A and wb behave like p" 1 at high densities, 
and thus C(w a ,wb) — 1 to first order. We then find 



a 2 (9)w A (9)w B (9)p(9)d 2 e 



(54) 



This expression should be compared with Eq. (|3*6l . 

If instead p — > + , then <3(sa 5 Sfl) — > 0~, Y(s A , sb) — > 
1 and thus C(iOyi, wb) — > v/{w a wb)- Note that here we 
are assuming ^ and wb{9) ^ 0. In the same 

limit, we have [see Eqs. (|19|) and (|3L)fl ] 



p(6/)d 2 6» . 



Hence, we finally find 



cr 2 {9)p{9)d 2 9 



p(9) d 2 9 



(55) 



(56) 



This result is consistent with the upper limit for T a ob- 
tained in Sect. 13.231 

4. Conclusion 

In this paper we have studied the statistical properties of 
a smoothing technique widely used in Astronomy and in 
other physical sciences. In particular, we have provided 
simple analytical expressions to calculate the average and 



the two-point correlation function of the smoothed field 
f(9) defined in Eq. The results generalize what was 
already obtained in Paper I and II to the case where ob- 
servations are carried out in a finite field, with a non- 
uniform spatial density for the measurements, and with 
non-uniform measurement errors a 2 (9). These generaliza- 
tions together greatly widen the range of applicability of 
our results in the astronomical context. Finally, we have 
shown several interesting properties of the average map 
and of the two-point correlation function, and we have 
considered the behavior of these quantities in some rele- 
vant limiting cases. 
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Appendix A: Alternative derivation 

In this appendix, we derive the same results obtained in Sect. [5] 
using a more direct method. Although not necessary, this al- 
ternative derivation is helpful in order to fully understand the 
whole problem and also clarifies some of the peculiarities of 
the equations derived in Paper I (cf., in particular, the case of 
vanishing weights). 

The derivation will follow quite closely the one adopted in 
Paper I, with the needed modifications due to the finite-field 
and the non-constant density. The only significant exception 
will be the use of a different strategy in performing the so- 
called "continuous limit" (because of the finite field, we cannot 
perform the limit N — > oo, but we must rather take N as 
a random variable). Note that throughout this appendix we 
will drop everywhere the index A, so that, e.g., wa(0) will 
be written just as w(9). This simplification should not create 
ambiguities, since anyway here we are concerned only with the 
value of / at 9a- 

Let us consider a field £1 and locations randomly distributed 
on this field with density p(6). Let us assume, for simplicity, 
that w(0) is strictly positive in f2; if this is not the case, we can 
always redefine fi to include only points inside the support of 
w (cf. discussion in Sect. The expected average number 

of locations in is given by 



N 



P (o) d 2 e 



(A.l) 



The actual number of points will be a random variable follow- 
ing a Poisson distribution with average N. In reality, since we 
are accepting only cases where there is at least one point in- 
side the support of w (we could not define / otherwise), N will 
follow the probability 



PAT (JV) 



N 



1 



-n jv! 



(A.2) 



where N is assumed to be a positive integer. Note in particular 
that the normalization factor takes into account the "missing" 
N = probability. 

Since the locations are distributed inside fl according to the 
density p(9), a single location follows the probability distribu- 
tion p(0)/N; note that the factor 1/N is needed here in order 
to satisfy the normalization of probabilities (the integral on 
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Q must be unity). Hence, the probability of having exactly TV In order to further simplify the definition of C, we use a tech- 
locations inside at the positions {9 n } (with n £ {1, . . . , TV}) nique similar to the one adopted in Paper I. Namely, we define 



is given by 

N 

P({On}) = PN (N)H 



P(8n) 

TV 



1 TV! 



(A.3) 



We can now use this probability distribution to evaluate the 
expectation value of / defined in Eq. @: 



(/> = E 



d 2 e N p({0 n }) 



ELM***) 



i 00 i f f 



jv/(giMgi) 



Similarly to Paper I, we now define, for each TV, the random 
variables 



vn = E W(0 n ) 



(A.5) 



In order to evaluate the associated probability distribution, we 
consider the probability distribution for values of the weight w 
(i.e., we use Marcov's m ethod; see, e.g.. IChandrasekharlll943l : 
iDeeuchi fc Watsonlll987h . The probability density for a loca- 
tion to have a weight w, in particular, is given by 



Pw (w) = I d 2 e P (0)s( w -w(o)) 



(A.6) 



This allows us to evaluate the probability distribution for j/jv 
as 



TV 



i^r / d 2 e 2P (o 2 y-- [ 

Jn Jn 



X 5{y N - w(0 2 ) 



'(On)) 



dw 2 p w (w 2 )- ■ 
x S(y N - w 2 



dWN Pw(vjn) 



W N 



Using this probability, we can rewrite Eq. 1A.4II as 



</} = 



7£ 



TV! 



d 2 9! p(0i)w(0i)f(0i 



f 

Jo 



N=l 

Np VN (yN)dy N 



w(Oi) + yjv 

The form of this expression justifies the definition 



C(w) 



TV 



JV-l foo 



TV! 



N p VN (yN)dy N 
w + yN 



(A.8) 



(A.9) 



so that we have 

(/)= / d 2 8 1 p(G 1 )w(e 1 )f(e 1 )c{w(e 1 )) 

Jn 



W, the Laplace transform of p u 



W(s) = C[p w ](s) = 



"p w {wa) dw 



1 

iV 



-8W(0) 



p{0) d 2 e 



(A.ll) 



We similarly define, for each AT, Yn, the Laplace transform of 
Pun ■> ^ 



Y N (s) = C\p yN ](s) = 
= [W(s)] N - X 



Py N (VN)e ayN dy N 



(A.12) 



We now use the following property of Laplace transforms: if / 
is any function, and xo any positive real number, we have [see 
Eqs. (14), (20), and (25) of Paper I for a proof] 



(A.4) C[C 



/(*) 



da; 



Xo + X 

Using this in Eq. JA.9II we find 
1 ^ TV 



(A.13) 



cm 



T N-1 



JV=1 



TV! 



N£[y N ](w) 



-C 



E 



TV 



JV-l 



W 1 



(w) 



e N - 1 
1 

e^ ~ 1 
1 



-C 



N v 



(w) 



l-e-ff 
Finally, we define Q as 

Q(s) = NW(s) -N = 

so that we finally have 



L[e NW }(w) 



r r NW-N1 , \ 



i]p(e)d 2 e. 



C(w) 



tC c 



l-e- N 
This completes our proof. 



» 



(A.14) 



(A.15) 



(A.16) 
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